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Abstract 

Simulations on a Lennard-Jones computer glass are performed to study effects 
arising from defects in glasses at low temperatures. The numerical analysis 
reveals that already a low concentration of defects may dramatically change 
the low temperature properties by giving rise to extrinsic double-well poten- 
tials (DWP's). The main characteristics of these extrinsic DWP's are (i) 
high barrier heights, (ii) high probability that a defect is indeed connected 
with an extrinsic DWP, (iii) highly localized dynamics around this defect, 
and (iv) smaller deformation potential coupling to phonons. Designing an 
extension of the Standard Tunneling Model (STM) which parametrizes this 
picture and comparing with ultrasound experiments on the wet network glass 
0-B2O3 shows that effects of OH-impurities are accurately accounted for. This 
model is then applied to organic polymer glasses and proteins. It is suggested 
that side groups may act similarly like doped impurities inasmuch as extrinsic 
DWP's are induced, which possess a distribution of barriers peaked around a 
high barrier height. This compares with the structurlessly distributed barrier 
heights of the intrinsic DWP's, which are associated with the backbone dy- 
namics. It is shown that this picture is consistent with elastic measurements 
on polymers, and can explain anomalous nonlogarithmic line broadening re- 
cently observed in hole burning experiments in PMMA. 



I. INTRODUCTION 



Unlike simple systems such as atoms, nuclei or crystals, complex systems cannot be 
entirely characterized by the concept of a ground state and elementary excitations: The 
ground state of systems such as glasses, spin glasses or proteins is highly degenerate, and in 
place of energy levels or quasiparticles, we must speak in terms of a high- dimensional "en- 
ergy landscape" in configuration space.EFi Basically the dynamics of the system is mapped 
on transitions between different basins of the energy landscape. Due to its complexity, 
this landscape is very rough and only very qualitative features are known about it. Ex- 
periments on myoglobin can be interpreted by saying that slightly different structures — 
conformational substate (CS), separated by high barriers — are organized in hierarchical 
tiers, ill Glasslike behavior is also predicted for folded proteins (tertiary structure) by both 
theory and experiment. §00 Somewhat closer information about the multiple minimum en- 
ergy landscape also comes from computer simulations on proteins, atomic clusters, and 
structural glasses.00'0 

Whereas experiments at elevated temperatures (around the glass transition temperature 
Tg) are sensitive to the overall topology of the energy landscape, experiments far below Tg 
are sensitive to its fine structure. In the framework of the standard tunneling model (STM) 
this fine structure has been successfully characterized by using a statistical approach.lliini 
The STM explains many properties of glasses in the Kelvin regime. It assumes that via local 
rearrangements the system can switch between pairs of local energy minima. Formally, this 
corresponds to the dynamics of a particle in a double- well potential (DWP). At temperatures 
in the Kelvin regime the dynamics has to be described in quantum-mechanical terms. The 
resulting tunneling dynamics only involves the lowest two energy levels. Hence the DWP's 
act as two-level systems (TLS's). An important aspect of the STM is the assumption of 
a broad distribution of parameters characterizing DWP's and hence TLS's. In particular, 
this broad distribution implies that transition processes occur on a broad spectrum of time 
scales. 



1 



Most low temperature thermal, acoustic and optical properties of glasses can be explained 
in this picture.lll However, recent hole burning experiments on a polymer glass (PMMA) in 
Haarer's group report a systematic disagreement with the STM.lii^ The authors found a 
logarithmic time dependence with a crossover to an algebraic behavior after ca. 3 h. At 
the same time equivalent measurements have been performed on a protein (myoglobin) in 
Friedrich's group0 with the result of almost no hole broadening up to 3 h and an algebraic 
time dependence thereafter. Very recently these results have been extended to longer waiting 
times.@ Though the logarithmic behavior is in agreement with the STM, the algebraic 
behavior is not. 

Apparently unrelated anomalies can be observed in sound attenuation experiments. In 
undiluted network glasses (Si02, Ge02) a single relaxation peak appears around 50 K 
which is consistent with a broad distribution of DWP's.ilhii However, an additional di- 
lution with OH-impurities gives rise to the existence of a second relaxation peak at higher 
temperatures.0J2i Specific mechanisms have been proposed giving rise to DWP's connected 
with the OH-impurities.ii 

Interestingly, sound attenuation measurements on various polymers also show an addi- 
tional strong increase at temperatures above ~ 30 — 50 K.&^ Different absorption peaks — 
labelled by Greek letters (a,/?, 7,5) in appearance of decreasing temperature — have been 
observed in the temperature regime between ~ 30 K and room temperature.@'0 Usually, 
the low temperature peaks (5, 7,/5) are attributed to side chain motions, which perform 
thermally activated reorientation processes even far below Tg.&iHli3 

This similarity between sound absorption in OH-doped network glasses and in polymer 
glasses may justify to think of side chain as "defects" , because both hydroxcyl groups and 
side groups give rise to well-defined localized DWP's. Thus one may ask the question whether 
the statistical approach of the STM still holds as soon as a significant fraction of DWP's is 
related to the dynamics of some well-defined defects. 

The conjecture we want to investigate in this paper is: (1) the existence of minority 
components of different structure (impurities as compared to the network or side groups as 



compared to the main chain) give rise to additional locahzed dynamical degrees of freedom; 
(2) associated with these degrees of freedom are induced or extrinsic DWP's which possess 
a nonuniversal, i.e., different from the STM, distribution of parameters; (3) these extrinsic 
DWP's cause a second absorption peak of sound at higher temperatures and a nonlogarithmic 
hole broadening at ultralong time scales. In what follows both OH-impurities and side groups 
will be referred to as defects. 

Rather than just postulating the nonuniversal properties of defects, we first show via 
computer simulations of a simple model glass that the presence of defects indeed has dramatic 
consequences for the distribution of DWP's. We choose a binary Lennard- Jones (LJ) glass 
which has been tested in previous simulations and has proven to confirm the assumptions 
of the STM.@ The LJ parameters are chosen to represent NiP. To this LJ glass we add 
"defects" (minority component) in form of a third LJ particle. The size of the defect will 
be controlled via its LJ parameters by a variable interaction length. The main finding of 
these simulations is that, in case that the minority component is somewhat smaller than 
the original glass components, they give rise to additional or extrinsic DWP's with high 
energy barriers whereas in case of large defects no extrinsic DWP's occur. The existence 
of small particles also has some impact on the intrinsic DWP's which are connected with 
the majority component. Although these results might not be representative for all glasses, 
one can speculate that as a general result the distribution of DWP's in systems with defects 
can be viewed to additionally comprise typical high barriers. Based on this observation, we 
design a model which extends the STM, and compare to experiments on OH-diluted network 
glasses, polymer glasses and proteins. 

The outline of this paper is as follows: In Sec. II, we describe the simulation method 
and present the results; in Sec. Ill, an extension of the tunneling model is proposed taking 
into account an additional contribution of DWP's with high barrier heights; in Sec. IV, we 
compare with experiment, and, in Sec. V, we conclude. Mathematical details are relegated 
to the Appendix in order not to interrupt the main flow of the paper. 
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II. SIMULATIONS 

A. Model for the computer glass 

The model LJ-type glass has been adapted from the work of Weber and Stillingerll and 
is devised to represent the binary glass former NiP. The system contains 80% type 1 (Ni) 
and 20% type 2 (P) particles. The interaction is pairwise and reads 

Aij[{aijr / oY^"^ - 1] exp[-o-/(a - aj^r)], for < atjr <a 

Vij = \ (1) 

else 

with a = 1.652 cr where a = 2.2 A is the unit length. Values for Aij and aij are An = 8200 
K, an = 1.0 (Ni-Ni), A12 = l.5Au,ai2 = 1.05 (Ni-P), A22 = 0.5 An, and 022 = 1-13 (P-P). 
The mass density po is 8348 kg/m^. We have simulated systems with = 150 particles 
and use periodic boundary conditions. We have checked that within statistical errors a 
simulation with = 500 particles yields identical results. The equilibrium distance between 
two particles of type i and j is given by parameters 2^^^ a /aij. In extension to the above 
model we include a single particle which we formally denote type 3 particle. Its potential 
parameters are chosen as a^a = a and Ai^ = An. The value of a determines the equilibrium 
distance between this particle and the rest of the glass and hence characterizes the effective 
size of this particle. The simulations have been performed in the range 0.7 < a < 1.6. Of 
course, one expects that for a ~ 1 this particle basically behaves like the rest of the glass 
for which 96% of all pair potentials are characterized by a = 1 or a = 1.05. 



B. Simulation method 

In previous work Heuer and Silbey developed a systematic search algorithm for DWP's.iil 
Starting from some T = configuration of a glass the idea is to search the energy landscape 
around some local minimum for a second local minimum. This pair of minima can be 
viewed as a DWP. The full cooperativity of the transition between both minima is taken 
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into account. In general all particles will be involved in the transfer between both minima. 
Let di denote the transfer distance of the i-th particle. The total distance of both minima 
in configuration space can then be expressed as 

^^-E4 (2) 

i 

The algorithm takes advantages of the observation that the existence of DWP's is related 
to a somewhat localized motion of particles. For different initial T = configurations we 
analyzed the region around the defect and checked whether we can find a DWP close to this 
defect. This is achieved in three steps. In a first step we choose the defect particle as well 
as its 15 nearest neighbors. Then in a second step we minimize the potential energy via 
some simulated annealing routine by moving the 16 selected particles such that the distance 
in configuration space to the original local minimum is = OA a. By this condition we 
avoid that the system just moves back to the original minimum. Finally, we relax the whole 
system to find a new local minimum. In case that close to the original minimum there exists 
a second minimum the algorithm will likely find this new minimum. A careful analysis for 
this choice of d^ has shown that DWP's with do < d < a are systematically found by this 
algorithm. Anyhow, the fraction of DWP's with d < OA a is very small. In present work 
we restrict ourselves to DWP's with d in this range. For the binary LJ glass it turned out 
that there exist one DWP per 100 particles. Hence in a small region around the defect one 
expects at most a single DWP. Further details about this algorithm as well as a discussion 



of its limitations can be found in Ref. 32. 



We wish to distinguish the properties of extrinsic and intrinsic DWP's. For this purpose, 
we determine for all DWP's which particle moves most during the transition between both 
local energy minima and denote this distance by ds- If this particle is the defect we denote 
this DWP extrinsic, otherwise intrinsic. Finally the effective mass is defined via 

p = dydl>l. (3) 

The value of p is a measure for the number of particles which are basically involved in the 
formation of the DWP. Finally, the potential height of a DWP is found by an explicit search 



for the saddle point. The present analysis is based on the analysis of approximately 3000 
independent initial configurations. 

C. Results of the simulations 

First we determine the probability that a DWP exists close to the defect. In Fig. 1 the 
probability is shown that the algorithm, described in Sec. IIA, finds an extrinsic DWP. One 
can see that the number of extrinsic defects dramatically increases for large values of a. For 
the smallest particles analyzed in our simulation runs nearly 80% of all defects are connected 
with a DWP. Of course, this number has to be viewed as a lower limit since it is possible 
that our algorithm may miss some DWP's. However, this already large number shows that 
the algorithm works realibly. For a = 1 the defect is basically identically to the rest of the 
glass. Hence the probability at a = 1 has to be identical to the probability (normalized 
per particle) to find a DWP in the whole glass. The present result agrees with the already 
above-mentioned finding of Ref. ^ that on average one DWP exists per 100 particles. For 
us it was very surprising that already for a = 1.2 a significant number of extrinsic DWP's 
can be observed. Hence only small differences between defect and structural glass seem 
to strongly enhance the probability for the formation of a DWP. For a = 0.7 no extrinsic 
DWP has been found. A huge particle is locked in a cage formed by the surrounding small 
particles without a chance to escape. 

In Fig. 1 the probability to find intrinsic DWP's is displayed, too. Originally, the simu- 
lation yields the probability that there exists an intrinsic DWP's having a significant spatial 
overlap with the initially selected particle cluster; see the description of the algorithm above. 
For a = 1 this probability is larger than the probability for the formation of extrinsic DWP's 
although no difference between intrinsic and extrinsic DWP's should exist. The reason is 
trivially related to the asymmetry in the definition of intrinsic and extrinsic DWP's. For 
a = 1 the ratio of the number of intrinsic and extrinsic DWP's should be close to the 
number of initially selected particles around the defect (here: 15 particles). This unwanted 
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statistical effect can be removed by scaling the probability of formation of intrinsic DWP's 
such that for a = 1 this probability is equal to that of the formation of extrinsic DWP's. 
This allows a direct comparison between intrinsic and extrinsic DWP's. All data in Fig. 1 
for the intrinsic DWP's are scaled in this way. Interestingly, the number of intrinsic DWP's 
close to the defect strongly depends on the size of the defect. The number is largest for 
q; ~ 1, but strongly decreases for larger or smaller defects. We interpret this observation as 
follows: for a very different to unity the defect does not participate in the dynamics of the 
glass because it does not move at all (small a) or because it only moves by itself (large a). 
Hence from the viewpoint of the glass particles there exists an "alien element" which forms 
a barrier for the surrounding particles. Thus close to the defect the degrees of freedom and 
hence the chance to form a DWP are significantly reduced. Only for ct 1 this effect is 
irrelevant since the defect may simply participate in the dynamics. 

In order to characterize the nature of the extrinsic DWP's somewhat closer, we plot the 
average effective mass p of the DWP's in dependence of a; see Fig. 2. One observes a 
transition from cooperative to strictly localized dynamics. Already for a = 1.2 the effective 
mass of extrinsic DWP's starts to decrease. For a = 1.6 the extrinsic DWP's are related 
to motions of the defect alone. In the hmit of very small particles this result could have 
been expected. Basically the defect jumps between different interstitial positions formed by 
a fixed environment. It is likely that for large a a defect has more than one direction in 
order to find a second minimum. From the present simulations this suggestion cannot be 
quantified. 

In Fig. 3, we plot the distribution of barrier heights of intrinsic DWP's with the presence 
of a close-by defect as well as the distribution for extrinsic DWP's related to a = 1.6. In 
both cases a broad distribution of barrier heights is observed refiecting the statistical nature 
of the formation process of a DWP. However, obviously the extrinsic defects in case of small 
defects are significantly shifted to higher potential heights. This is quantified in Fig. 4, 
where we plot the average barrier height V in dependence of a. In agreement with Fig. 
3, one observes a strong increase of the value of the average potential height with a. As 
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already discussed in Ref. |32|, the ability to cooperative motion, i.e., a large effective mass 
p, tends to decrease the necessary energy changes along the way between both minima of a 
DWP. In the opposite limit p = 1 a. single particle has to move by itself, hence experiencing 
strong resistance by the environment leading to a larger value of the barrier height. This 
general statement is again confirmed by the observed correlation between the effective mass 
and the average barrier height; compare Figs. 2 and 4. 

Finally, in Fig. 5, we show the value of the deformation potential 7, normalized such 
that 7(0; = 1) = 1. The theoretical basis for this calculation can be found in Ref. |3^. 
The deformation potential decreases with increasing a. This effect can easily be explained. 
The number of nearest neighbors of the defect strongly depends on its size and decreases 
with decreasing size. This can be easily visualized for a regular hexagonal lattice in two 
dimensions where particles on interstitial sites are surrounded by three particles whereas 
the coordination number of a regular lattice site is six. Since the deformation potential is 
related to the interaction of the defect with the environment it is not surprising that the 
deformation potential strongly depends on the number of nearest neighbors. 



III. EXTENSION OF THE TUNNELING MODEL 

In the STM it is assumed that the asymmetry energy between two wells, e, and the 
tunneling parameter, A = {d/h)y/2mV , are uniformly distributed over a wide range, 

P{\,e)d\de = PodXde (4) 

with a constant Pq ~ 10"^^ — 10'*^ J^^ m^^ for most glasses. Here, m represents the mass 
of the tunneling unit, and d the distance between the minima of the double-well potential. 
The tunneling parameter is related to the tunneling amplitude. A, through a WKB relation 

A = huJoe~^, (5) 

where huo is essentially the zero-point energy, and to the one-phonon relaxation rate by 



8 



where 

A^in(^) = \og{huo^/E). (7) 



The detailed form of the maximum rate, i?max at fixed TLS-energy, E = yAM-?, 
depends on the spectral density of vibrational modes in the energy landscape. Usually, a 
deformation potential coupling, with parameter 7, to a Debye-spectral density, J{uj) oc uj^, 
is imposed. The one-phonon rate then reads R = {A/ E)'^J{E /h) coth{E/2 ksT) which 
together with (j^) and (||) provides 

^max(^) = cT^a;^ cothx , (8) 

where x = E/2kBT and c = 7^ {2kB)y27ih^gv^. Here, g is the mass density and v the 
averaged sound velocity. The ensuing hyperpolic distribution for A and R, P{e, A) = 
Po/A or {e/E)P{E,R) = Pq/2R, explains in particular the broad distribution of relaxation 
processes in glasses on a logarithmic time scale. 

The simulations suggest to extend the tunneling model by adding a second peak for a 
glass with additional defects (OH-groups) or a polymer glass (side chain motion) [cf. Fig. 1]. 
In order to restrict ourselves to a minimum set of new parameters, we only characterize this 
second peak by its mean value and its width in the following model: 

P(e,A) = Poe-(^-^°)'/2ag ^ g-(A-A0V2a? ^ ^ > x^^(^e), (9) 

where Pq and Pi weights the density of the intrinsic and extrinsic TLS's. By formulating 
the new distribution for the tunneling parameter. A, we additionally assume some relation 
between the potential height V and the distance between the two minima of the DWP, d. 
The standard assumption \s,V (x d?. For Pi = and o"o ^ 1 the standard model emerges.00 
The parameter Aq determines the most probable tunneling matrix element. For all practical 
applications this parameter in irrelevant and will be set to zero hereafter. A model including 
only the first term has previously been used by Jankowiak and Smalhil As we will see below, 



a combination of both terms is needed in order to explain several experiments consistently. 
These authors have also emphasized the importance of extrinsic TLS, which they attribute 
to the presence of the chromophores in the probe. 

The ensuing distribution in A = hujo then reads after neglecting terms depending on 
(l/2CT|)log(Ai^„/A)«l 



where P{ = Pie and 



KA) = ^ - ^ log(/ia;o/A). (11) 
(Tf /erf 



Identifying ctq with the parameter Amax of the STM, one deduces that 10 ~ ctq ~ 30. The 
value of Pi depends on the concentration of the defects. For a typical concentration of 
OH-groups in network glasses, one expects Pi ~ Pq. In polymer glasses, however, the 
number of extrinsic DWP's can be larger because of the large number of side chains. If 
every side chain contributes one extrinsic DWP, and if only one out of ~ 100 monomers of 
the backbone forms an intrinsic Dwpii, we find as an upper bound that Pi ~ 100 Pq. More 
difficult is the estimation of mean value and the width of the second peak. For a = 1.6 one 
observed numerically an average value of d ^ 0.6 a. Together with the average potential 
height Vav = 800 K and using the average atomic mass of NiP one may estimate Ai ~ 35 
and cTi ~ 25. However, for any realistic applications Ai (or Vav) and ai depend on the 
exact nature of the defects and hence should be treated as adjustable parameters. Due to 
the observation of several absorption maxima below the glass transition temperature Jh^IhI 
a more realistic model should contain more than one additional peak in the distribution 
of barrier heights. Since we are interested in low temperature properties, it is sufficient to 
include only the lowerst addition barrier peak in our model. 



IV. COMPARISON TO EXPERIMENT 
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A. Sound attenuation in OH-doped network glasses and polymers 



A possible method to measure sound absorption is to clamp a small glassy plate on one 
end and to drive it electrostatically to forced vibrations on its free end, the so-called vibrating 
reed technique.0^0 The measurement of the resonance frequency and the amplitude of the 
plate then determines the acoustic loss and the variation of the sound velocity. 

Measurements of the elastic properties of glasses — for example the internal friction, 
— have proven to be a powerful method for the study of tunneling states below ~ 5 K, 
where it is assumed that relaxation occurs with the one-phonon rate (P). At higher tem- 
peratures, typically between 10 K and room temperature, the elastic measurements pro- 
vide information about the distribution of potential barrier heights for thermally activated 
reorient at ions ,0 

with To ~ 10"^'^ s, of atoms and molecules. Assuming a simple oscillator model for the 
individual wells of a DWP, we have V (x dP and 

A = ^e-\ (13) 

TT 

where Eq is the zero-point energy of the oscillator, and 
With this, the distribution can be rewritten yielding 

P(e, v) = ^ e'^y~y^fi^~^l + ?!. e-'y-y^f/^-^i (15) 

with Vi = XiEo and 5, = aiEo, and the internal friction, Q^^, can be calculated with the 
formalism developed in Ref. BTl. 



Comparing the data for dry and wet a-B203 in Figs. 7 and 8 of Ref. one clearly 
sees that the sound absorption is strongly effected by the presence of OH- impurities. For 
the dry sample a broad absorption peak is observed around ~ 50 K, while for the wet 
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sample a shoulder between 30 and 80 K is followed by a narrow peak at ~ 160 K. Similar 



observations have been made in Ref. where the dry a-B203 sample provided one broad 
peak whereas the wet sample showed two maxima. The authors in these references concluded 
that intrinsic DWP's cause a loss peak at lower and OH-induced, i.e., extrinsic DWP's a 
peak at higher temperatures. Independent of any theoretical framework, these experiments 
directly indicates that the barrier heights of these OH-related DWP's must be higher than 
those of DWP's related to the glass itself. Furthermore, the mere fact that approximately 
1% OH-content can have such a significant influence on the relaxation properties shows that 
OH-defects have to be very efficient in prompting the formation of DWP's. 

In Fig. 6, we have plotted the internal friction data of wet 0-6203 sample together with 
a fit from our model (|l^) andU 

= fdefdV P(e, V) sech\e/2kBT) (16) 

where r is given by the one-phonon rate or the Arrhenius rate (|T2|) depending on whether 
the temperature is in/below the plateau regime or above, respectively. The data are from 
Ref. The solid line shows the contribution of intrinsic and extrinsic DWP's; the dashed 
line the contribution of intrinsic DWP's. The parameter values are given in the figure 
caption. We see that elastic measurements are a powerful tool to determine the barrier 
heights distribution, and that our computer simulation based model (^) accounts well for 
the infiuence of OH-defects in network glasses. 

Let us now discuss the elastic properties of polymers. 

Nittke et al.S have recently measured the low temperature (< 50 K) elastic prop- 
erties of the polymer glasses polymethylmethacrylate (PMMA) and polystyrene (PS) at 
uj = 2n X 240, 535 and 3200 Hz, following earlier measurements by Federle and Hunklinge 
at = 27r X 15 MHz and Geis, Kasper, and Hunklinger^ at u; = 27r x 430 Hz in PMMA. 
In the latter work, also the high temperature absorption between 50 and 300 K has been 
measured following previous work, cf. Refs. pB|-p5|. Below ~ 30 K the acoustic loss of 
PMMA and PS shows the typical temperature dependence of network glasses like Si02: the 
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plateau merges in the thermally actived relaxation peak around ~ 5 K. The missing first 
relaxation peak of the acoustic loss in PMMA at uj = 2it x 535 Hz (cf. Fig. 12 in Ref. 



30| or Fig. 4 in Ref. ^) can well be explained by the low applied frequency uj/2tt and a 
small value for the cutoff, \4iax, in the barrier distribution of the STM-states.S However, 
above ~ 30 K both polymers behave anomalously compared to (undoped) network glasses 
by exhibiting a strong increase in their acoustic loss. Early measurements show that this 
increase is characteristic for polymers and continues up to the glass transition temperature 
with various (5, 7,/3) relaxation peaks superimposed to it.iiill Recently, multidimensional 
NMR measurements have proven that the molecular origin of the (3 relaxation in PMMA 
at 330 K (10 Hz) is a large- amplitude flip of the methacrylate side group around the C-C 
bond and a concomitant main chain torsion.ill From the evidence at present available the 
nature of the groups responsible for the low temperature 6, 7 peaks which have been ob- 
served in different polymers is less clear, however, they have traditionally been attributed 
to thermally activated rotation of sub groups within the side chains. Examples are the n- 
alkyloxycarbonyl side group rotations of the polymethacrylates, see Figs. 3 and 4 in Ref. 
In PMMA these side groups do not exist and absorption peaks observed at temperatures 
below the P peakil are thought to be due to the presence of small amounts of dibutylphta- 
late (DBP) — a plasticizer — or higher polymethylacrylate esters. Studies on purified samples, 
which revealed a pure exponetial increase, Q^^ oc e""^, between ~ 30 and ~ 300 K, support 
this interpretation.il0 However, even with a broad distribution, it is not clear whether the 
exponential increase is connected with the (3 process. It is conceivable that further molecular 
dynamics occurs in the intermediate temperature regime. 

In any case we learn from these experiments that (1) the similarity between sound ab- 
sorption in OH-doped network glasses and in polymer glasses may justify to think of side 
chain as "defects" , because both hydroxcyl groups and side groups give rise to well-defined 
localized motion within DWP's with high activation barriers, and thus constitute additional 
degrees of freedom; (2) the energy landscape of polymers is structured at higher barrier 
heights: in addition to the uniformly distributed low barriers being responsible for the 
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plateau region below ~ 6 K, several tiers of high barriers with typical activation energies 
can exist depending on the chemical compound of the polymer. 

B. Hole burning in polymer glasses and proteins on ultralong time scales 

In a hole burning experiment the probe is doped with a large number of guest molecules 
having resonance frequencies in the optical range. Due to the disordered structure of the host 
each chromophore experiences a slightly different environment which detunes its resonance 
frequency; the absorption line becomes inhomogeneously broadened, accordingly. A fraction 
of these chromophores is photochemically or photophysically removed from the absorption 
line by laser irradiation. The inhomogeneous line then shows a narrow hole at the laser 
irradiation frequency. Each chromophore is interacting with local rearrangements of atoms 
in the host which conveniently can be described in the TLS picture. As a result, the 
resonance frequency of the guest molecule performs spectral jumps and its linewidth becomes 
broadenend. This process, which is called "spectral diffusion" , results in a partial refillment 
and, thereby, in a broadening of the hole. Hence, the time dependence of the hole width is 
a mirror of the relaxation processes in the host, and, accordingly, a tool to investigate the 
local dynamics of the host. In proteins and glasses the results can be interpreted by the 
concept of an energy landscape. 

In recent years this technique has been used extensively by various groups on glasses 
and proteins.i'@^'&^l Most measurements on glasses for times t < 3 h confirm the STM, 
which predicts according to {e/ E)P{E, R) = Pq/2R and 



of the hole width after some reference time to.E3'E3E^ Here, (C) is the chromophore-TLS cou- 
pling strength. This accords with the broad distribution of relaxation times on logarithmic 




(17) 



a logarithmic growth 




(18) 
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time scales. 

Recently, Haarer and coworkeJlUll have performed hole burning on purified PMMA 
at 0.5, 1 and 2 K up to extremely long times (from 10 s to 10 days), see data points in 
Fig. 7. They found a log(t)-dependence up to ca. 3 h, but a v^-behavior between 3 h and 
10 days. Furthermore, Friedrich and coworkerlii0 have performed equivalent measurements 



on proteins at 100 mK and 4 K. In Ref. |I9|, these authors did photochemical hole burning 
on a glycerol/ water glass over 44 h with and without the protein myoglobin doped into the 
sample, see data points in Fig. 8. Recently, they also studied the protein cytochrome c in 
a glycerol/dimethylformamide glassS for 300 h, see data points in Fig. 9. The result was: 
(i) spectral diffusion broadening is always less in the protein than in the glass; (ii) the glass 
showed a logarithmic time dependence within the first 3 h, while there was practically no 
relaxation in the protein over this period; (iii) after about 3 h there is a strong nonlogarithmic 
increase in the line broadening for both the glass and the protein; (iv) contrary to the glass, 
the protein showed no aging effect indicating a gap in the relaxation rate between 1 s and 
~ 3 h. The protein and the glass data could be fitted with an empirical ad hoc ansatz for 
the distribution function 



P(e,A)=Po 



1 A 
A ^ A2 



A = constant, (19) 



resulting in a growth composed of a superposition of a log(t)- and a -\/t-term. For the 0.5 
and 1 K data in PMMA a single consistent parameter set with A/ks = 10^^ K has been 
used.0'lll For the 2 K data in PMMA, however, a 60% smaller value for the weight parameter 
A was needed.El For the protein data the question rises whether there is a log()f:)-term or 
not as the data can already be reproduced by the y^-term. This observation motivated 
the authors in Ref. ^ to discuss a possible alternative interpretation of the nonlogarithmic 
dynamics which is related to classical diffusion in the energy landscape. 

In both experiments the distribution function (|I9|) has been motivated by recent pub- 
lications focusing on the interaction of TLS's in glasses,&^ in particular by the scenario 
of coherently coupled pairs of TLS's, which was invented by Burin and Kagan.ill In this 
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formulation the second term in (|T9]) stems from pairs of single TLS's which are coupled by 
a (resonant) up-down transition. If the single TLS's are distributed as in the STM [first 
term in Eq. ([T9|)], a 1/A^-distribution emerges for the pairs. In Ref. |5^, a quantitative test of 
this picture has been performed with the conclusion that this strongly modified TLS model 
could, in principle, explain the deviation from the log(t)-dynamics. At higher temperatures 
pairs are expected to break up due to thermal fluctuations which destroy the coherence of 
the coupling. It has been suggested that this is the physical reason for the decrease of the 
weight factor A seen in the PMMA data at 2 K.El 

In the following, we will propose an alternative explanation, which is motivated by the 
acoustic measurement on polymers, and has no restriction towards higher temperatures. 

Instead of considering TLS-fiip-fiop processes as the physical cause of spectral diffusion, 
we take barrier crossings events. On the time scale of seconds and larger, both picture are of 
course equivalent due to the big ratio e/A of the relevant TLS or DWP's. We assume that at 
1 K the barrier crossing process occurs via tunneling with the one-phonon rate, Eqs. , 
and consider the tunneling parameter A, which in the simple oscillator model (|T^) is equiv- 
alent to the barrier height V . Because of the observation that the broadening significantly 
increases at 1 K between tc = 10^ and 10^ s, and that the distribution -P(A) oc 1/A^ well 
accounts for this phenomenon, we learn that the distribution P{X) increase exponentially 
around some value Ac ^ 1. We may estimate this value from the relation R{lK)tc ~ 1 
which provides with i?max(lK) ~ 10^° s^^ that Ac — Amin(lK) ^ 15. A natural way to 
interprete this increase in A is to assume that the energy landscape of a polymer glass is 
not structureless, but that it comprises high barriers in addition to structurelessly distributed 
lower harriers within each basin — as observed in the acoustic experiments on polymer glasses. 
This is the picture which Frauenfelder suggested for the energy landscape of proteins.&i'i We 
may identify transitions between deep basins with degrees of freedom (rotations or switch- 
ing) located at the side chains of the polymer glass or protein. The STM-like states may be 
associated with the more collective dynamics along the backbone or the network. One might 
think that already the backbone dynamics of polymers gives rise to dynamics with high bar- 
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rier heights, since this dynamics is related to conformational transitions. However, it has 
been shown for some polymers that conformational transitions occur on the same time scale 
as the a relaxation and thus can be neglected at temperatures far below Tg.S The reason 
is that conformational transitions involve highly cooperative dynamics in contrast to more 
localized processes we are here interested in and which are still active below Tg. 

The difference between protein and glass then simply lies in the amount of disorder and 
organization in the energy landscape. Though deep basins separated by saddle points exist 
in both materials, in proteins the distribution of small (STM-like) barriers in the basins is 
restricted to smaller values, i.e., o"o (protein) ^ ctq (glass), or, more likely, itself composed 
of many well-separated distributions yielding the hierarchical picture of Frauenfelder with 
STM-like states on the lowest tier. Indeed, recent stimulated echo experiments by Thorn 
Leeson and Wiersma can be interpreted in this way.0 Hence, whereas many different tiers 
can exist in proteins, the glass is less organized: the structurelessly distributed STM-states 
extend to much higher barriers, and only the passing from one basin to the next gives rise to 
an (in the sense of the STM) anomalous behavoir. Hence, we suggest a connection between 
the low temperature (~ 100 K) relaxation peaks observed in elastic measurements with the 
nonlogarithmic line broadening observed in spectral hole burning after 10'^ — 10^ s in the 
temperature regime around 1 K. A crossing of barriers as high as ~ 1000 K by tunneling 
at 1 K is only possible because of the extreme long waiting time of 10'^ — 10^ s. At higher 
temperatures (~ 5 K) the one-phonon rate should gradually be replaced by the Arrhenius 



rate (p!2D, which will change the temperature dependence, but not the time dependence of 
the hole broadening. 

The simplest way to quantify our ideas is to calculate the hole broadening from Eq. (Bl) 



and compare with the experimental data of Refs. p!6Hlq and p!9| , p0| . In the Appendix we give 
details of the derivation. We find for typical values of the shortest time a hole can be read, 
to ~ 1 s, 

Ar(t) = ^(C) fcsT [Poaov^/^'H^) + Piaiv^/^^Hi^) ], (20) 
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where 



= erf (^^^^^^^) - ("^^' '^jf "°') . (21) 

Here, erf(a;) (erfc(x)) is the (complementary) error function, and K has been defined in Eq. 
( |A.6| ). The first term yields with erf (2;) {2/ ^/t:)z in leading order in (I/ctq) the standard 
result ([ll). If the barrier crossing is thermally activated with rate 
easily finds that within the simple oscillator model (|l^ and (^) 

(23) 



/«(t) = erff^^^MZ!£)') - ,,i(kBT\og{t,/r,r 



f^^\t) ^ erfc {^''XT'^^^) - (^^^^Cl^) ' (^'^ 
with cTj = -EoCj, which, as mentioned above, changes the temperature but not the time 
dependence of the line width. 

In Figs. 7 to 9, we compare Eq. (pOD - (|22| ) with the experimental data for PMMA, 
glycerol/water glass and myoglobin, and cytochrome, respectively. The parameter values 
are given in Table I and II. In the glasses PMMA and glycerol/water glass, logarithmic line 
broadening occurs below ~ 10"^ s, as predicted by the STM. The experimentally detected 
algebraic -\/t-time dependence applies approximately between 10^ and 10^ s. As is clearly 
demonstrated in Fig. 7, it arises from a superposition of the two terms in Eq. (pO|). Hence, 
the experimentally found exponent 1/2 has no physical significance in this model. The 
protein data can be fitted without a log(t)-term (i.e., Pq = 0) as the dashed lines in Fig. 8 
and 9 illustrate. In any case the contribution of the log(t)-term is very small as can be 
inferred from the dash-dotted curve in Fig. 9. This is consistent with the absence of aging 
observed in Ref. ^ (see data points corresponding to different equilibration times in Fig. 9), 
which indicates a gap in the distribution of relaxation times between ~ 1 s and 3 h. The 
hole broadening of glycerol/water glass and myoglobin embedded in the glycerol/water glass 
could not be fitted with the same set of parameters, which clearly indicates the protein 



comprises itself TLS's which dephase the chromophore. The numerical values for K, Po{C) 
and ctq are in the range known from other experiments. The numerical values for the new 
parameters Ai, ai, and Pi/Po lie well within the range we estimated from the numerical 
simulations. Using the parameter set for PMMA in Table I for the acoustic measurement of 
Refs. pP|J5(]| is difficult because we do not know the exact relation between V and A. However, 
assuming £'o(PMMA) ~ 10—20 K we see that the increase in the sound absorption in PMMA 
should be at a lower temperature than for the OH-doped Boroxid glass — as observed in the 
experiment. 



C. Comparison of elastic and optical properties of PMMA and experimental test 

Our picture of a structured barrier heights distribution seems to be consistent with elas- 
tic and optical measurements on polymers. However, in PMMA the puzzle remains that 
purified samples contain no sound absorption peak in the temperature regime below the (3 
peak (~ 330 K). Hence, no peaks in the barrier heights distribution seem to exist in the 
regime below say 1000 K, which is the relevant order of magnitude to obtain a nonlogarith- 
mic line broadening between 10^ and 10^ s. However, one must exercise caution since the 
experimentally observed exponential increase in the acoustic loss between 20 and 300 K can 
hide different dynamic molecular processes. Furthermore, in principle, our simulations show 
that any impurity can generate such an additional peak in the barrier distribution, and that 
defects are very effective in prompting external DWP's. A very low concentration of defects 
is sufficient to have dramatic effects. Hence, it is possible that the presence of impurities 
such as exchange gas atoms, monomers, fragments, hydrocyl groups, etc., is reponsible for 
this process. Also, it is unclear what the role of the chromophores may be, although the 
observation of a logarithmic line broadening in proteins, where the chromophore only re- 
places a group of a similar structure (i.e., does not change the local configuration of the 
protein), makes the presence of the chromophore a doubtful cause of the nonlogarithmic line 
broadening. Furthermore, our simulation did not show any additional peaks in the barrier 
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heights distribution function upon adding large impurities to the LJ-glass. 

There are several possibihties for an experimental test whether the physical cause of 
the nonlogarithmic line broadening is due to the explanation given in this work or due 
to the TLS-TLS interaction as pointed out in Ref. |5^. First, it would be interesting to 
know whether this effect is universal, i.e., whether hole burning on different polymers give 
identical results: same power law and same crossover time from the logarithmic to the 
nonlogarithmic (or algebraic) behavior. If this were true, it would hint towards TLS-TLS 
interaction as the physical cause. Even in this case, however, it would be difficult to argue 
in favor for this explanation at temperatures above ~ 5 K. It seems so that independently 
of a possible coherent TLS coupling around 1 K, the explanation given in this work is 
more intuitive at higher temperatures. Especially, the observation of a transition from 
(pip and (p2D to (p3D and p^) would favor this interpretation. Next, elastic measurement 
up to higher temperatures on exactly the same sample which has been used for the hole 
burning experiment, i.e., including the chromophore, would also give additional information. 
Finally, the question how polymer glasses without side chains, like polyethylene (PE) or 
polybutadiene (PB), behave in hole burning at ultralong time scales seems to be interesting. 
Following the interpretation of the present work, one would expect no or a later deviation 
from a logarithmic hole broadening (if ctq is sufficiently large) as compared to polymers like 
PMMA because of the lack of "defects." 



V. SUMMARY AND CONCLUSIONS 

We have analyzed the energy landscape of glasses containing defects. Specifically, we 
refer to OH-defects. However, the experimental results on polymers allow to speculate that 
also some side groups in polymers like the carboxyl group in PMMA may be viewed as a 
(generalized) defect. 

The analysis of LJ glasses has revealed that already small differences between defect and 
glass dramatically change the energy landscape around this defect, giving rise to additional 
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extrinsic DWP's. The main characterstics of these extrinsic DWP's are (i) high barrier 
heights, (ii) high probabihty that a defect is indeed connected with an extrinsic DWP, (iii) 
highly locahzed dynamics around this defect. As aheady discussed in the Introduction 
(i) and (ii) are in direct agreement with experimental observations on network-formers. 
Although the results are derived for a specific model glass we believe that this behaviour 
may be generic for other types of defect/glass-systems. 

This observation allows a straightforward extension of the STM which successfully ac- 
counts for the nonlogarithmic dynamics of hole burning experiments on polymer glasses 
and proteins. In terms of the energy landscape these low temperature experiments mainly 
probe their fine structure. Comparing the hole burning experiments for polymer glasses 
and proteins for short times leads to the conclusion that the fine structure is more distinct 
for polymer glasses. However, the presence of some high energy saddles, corresponding to 
localized motion, seems to be present in both systems. 
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APPENDIX A: 

In this Appendix we derive Eq. (^). Starting out from Eq. ([T7D, we first note that 

[e{E, R)/E]P{E, R) dEdR = [e{E, X)/E]P{E, A) dEdX 

= (Poe-("-"°)'/''^« + p,e-("-"^)'/^'^?) dEdX. (A.l) 

Rewriting Eq. (P) as 
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tR{E,X) = e-2[^-^mm(£;)-(l/2)log(ti?n.ax(E))]^ (^^2) 

we may use the step function approximation for the term (^(>-^oR{e,\) _ g-tR{E,\)^ -g^^ 

/ dAfe-*°^(^'^) -e-*^(^'^)) ... ^ / dA... , A.3 

where 

e(^) = max[A^i„(^), A^i„(E) + (1/2) log(toi?max(^))]. (A.4) 
For a hole burning experiment, we have to > 1 s such that at T ~ 0(1) K 

i{T) = A^i„(T) + (1/2) log(toi?max(r)) = (1/2) log(to^(T)), (A.5) 

where 

R{E) = i?„ax(^)e^^™^^^ = KTx coth x, (A.6) 

K = c{hujQ/2kBy [cf. Eq. (||)]. Note that, due to Eqs. (0) and (||), R depends hnearly on 
temperature. From the first term in Eq. ( [A.lD , we then find with0 



1 -(--)Wrf^=i(l + erf(i^)), (A.7) 



where erf(z) is the error function, that (after setting Aq = 0) 



POO 

AT^^\t) = (7rV3/i)(C)Po / dE sedi^{E /2kBT) 

Jo 



X 



From the second term in Eq. ([A.l|) , we find with ( |A.71 ) and the definition of the comple- 
mentary error function erfc(2;) = 1 + erf (—2;) that 



AT^^\t) = {ny3h){C)Pi / dE sec\i^ {E /2k bT) 

Jo 



X 



Upon noting that R{E) is approximately independent of E ioi E < 2kBT, we find Eq. (^ 
in the hole burning time regime (|A.5|) . 
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Figure captions 

FIG. 1. The probability to find extrinsic and intrinsic DWP's by our search algorithm. As 
discussed in the text the probabihty for intrinsic DWP's has been scaled such that for a = 1 
both probabilities agree. The observation that no extrinsic DWP's have been found for 
a — 0.7 is indicated by the arrow. Note the logarithmic scale for the probability. 

FIG. 2. The average effective mass p in dependence of a. 
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FIG. 3. The distribution of potential heights for a = 1.0 and a = 1.6. 



FIG. 4. The average potential height V in dependence of a. 



FIG. 5. The average deformation potential 7 in dependence of a. The data have been 
normalized such that 7(0 = 1) = 1. 



FIG. 6. Internal friction, Q ^, of the OH-doped network glass 0-6203. The data are from 



Ref. ^ The solid line is a fit with the distribution @; the dashed line line shows the 
contribution of the internal DWP's. The parameters are C = Po'y'^/gv'^ = 3.6 x 10"'', 
Eo/kB = 18 K, ao = 450 K, Pq/P^ = 0.22, Vi = 2000 K, and ai = 900 K. 



FIG. 7. Hole broadening in the polymer glass PMMA at 2 K (upper curve), 1 K (middle 
curve), and 0.5 K (lower curve). The experimental data are from Ref. |13,|15. The solid lines 
are fits with Eq. ( pUD with parameters given in Table I; the dashed line depicts the contribu- 
tion of the first term, the dash-dotted line the contribution of the second term in Eq. (pUD. 
A nearly algebraic line broadening between 10^ — 10^ min emerges from the superposition 
of both terms. 



FIG. 8. Hole broadening in glycerol/water glass (upper curve) and the protein myoglobin 
(lower curve). The experimental data are from Ref. |19[ The solid and dashed line are fits 
with Eq. ( pOD with parameters given in Table I and II. The protein data have been fitted 
without a log(t)-term (Pq = 0). 



FIG. 9. Hole broadening in the protein cytochrome c. The experimental data are from Ref. 



20| ; time between cooling and burning: 94 min (+), 4398 min (o), and 10093 min (*). The 
full and dashed line are fits with Eq. (|2^) with both terms (full line) and only the second 
term (dashed line). Parameters given in Table I and II. The dash-dotted line depicts a 
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possible contribution of a log(t)-term. 



Tables 



Table I: Parameters for PMMA, glycerol/water glass and cytochrome 

(intrinsic + extrinsic DAAT's): 





Ai 


o-Q 




Pi/Po 


Po{C) 


K [K-is-i] 


PMMA 


20 


30 




8 


0.95 X 10-^ 


5 X 10^^ 


glycerol/water 


20 


20 


n/Tj 


28 


3.49 X 10-5 


3.55 X 10^^ 


cytochrome 


19.7 


20 


1 


50 


0.35 X 10-5 


4.1 X 10^° 



Table II: Parameters for myoglobin and cytochrome 



(extrinsic DWP's): 





Ai 


^1 -Po Pi{C) 


K [K-^s-i] 


myoglobin 


20 


\/3 62.4 X 10-5 


6.6 X 10^° 


cytochrome 


19.7 


1 17.3 X 10-5 


4.1 X 10^° 
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FIG. 3: Heuer et al. 
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FIG. 5: Heuer et al. 
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FIG. 6: Heuer et al. 
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